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Ways to constrain neutron star equation of state models 
using relativistic disc lines 
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ABSTRACT 

Relativistic spectral lines from the accretion disc of a neutron star low-mass X-ray 
binary can be modelled to infer the disc inner edge radius. A small value of this 
radius tentatively implies that the disc terminates either at the neutron star hard 
surface, or at the innermost stable circular orbit (ISCO). Therefore an inferred disc 
inner edge radius either provides the stellar radius, or can directly constrain stellar 
equation of state (EoS) models using the theoretically computed ISCO radius for 
the spacetime of a rapidly spinning neutron star. However, this procedure requires 
numerical computation of stellar and ISCO radii for various EoS models and neutron 
star configurations using an appropriate rapidly spinning stellar spacetime. We have 
fully general relativistically calculated about 16000 stable neutron star structures to 
explore and establish the above mentioned procedure, and to show that the Kerr 
spacetime is inadequate for this purpose. Our work systematically studies the methods 
to constrain EoS models using relativistic disc lines, and will motivate future X-ray 
astronomy instruments. 

Key words: accretion, accretion discs — equation of state — methods: numerical 
— relativity — stars: neutron — X-rays: binaries 



1 INTRODUCTION 

The core density of a neutron star is typically 5—10 
times higher than the nuclear density. Such a dense mat- 
ter at a relatively low temperature (e.g., ~ 10 s K) can- 
not be probed by heavy-nuclei collis ion experiments or with 
observ ations of the early universe l|Lattimer and Prakash I 
120071 ) and references therein). Plausibly the only way to 
probe this degenerate matter is to constrain the theoreti- 
cally prop osed equation of state (EoS) models of neutron 
star cores (|Shapiro and Teukolskv"1ll983l ). For a given EoS 
model and an assumed central density, the stable struc- 
ture of a non-spinning neutron star can be computed by 
solving the Tolman-Oppenhei mer-Volkoff (TOV) equation 
jShapiro and Teukolskvlll983h . These stable stars trace a 
single curve in the mass (M) — equatorial-radius (R) space. 
Therefore, measurement of mass and radius of the same non- 
spinning neutron star would constrain the EoS models. If the 
neutron star spins, then the EoS models can be constrained 
by the reliable measurements of three independent parame- 
ters, such as mass, radius and spi n-frequency (t4 P i n ) of t he 
same neutron star (e.g., Fig. 1 of Bhattac harvva I (|2010h ). 
This is extremely difficult because of a number of unknown 
systematics. Low-mass X-ray binaries (LMXBs) can be par- 
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ticularly promising systems for such measurements, because 
several complem entary methods are av ailable for them. A 
recent review by iBhattacharvval (2010) describes how the 
simultaneous application of these methods can reduce the 
systematic uncertainties. 

One such method involves broad relativistic spectral 
lines from the inner portion of the accretion disc. The 
strongest among these fluorescent spectral emission lines 
is the one for the n = 2 — > n = 1 transition of the 
iron atom (or ion), and is observed from many accreting 
supermassive and stellar-mass black hol e systems (see 
iRevnolds and Nowak I (l2003i):[Mi!ler1 d2007l); and ref e rences 



therein). Recently. iBhattacharvva and Strohmaver I (|2007l ) 
has, for the first time, established that the broad iron 
lines from neutron star LMXBs also originate from the 
inner accretion discs. This discovery was soon confirmed by 
Cackett et al. I (|2008l ) using data from a different satellite. 
After these initial reports, the inner disc origin of broad 
iron line has been c onfirmed for se veral ot her neutron 
star L MXBs (e.g., iPandel et al~l (120081): ID Ai et al I 
d2D09h : ICackett et al. I J2009I); iPapitto et al. I J2009I); 



iReis et aL 

1 20091 ) ; ICackett et al. 



20091): Idi Salvo et al. I (|2009l ): llaria et al. I 



(20l3)) 



The broad iron line is affected by the strong gravity 
and spin of the compact star (neutron star or black hole). 
Before discussing how this relativistic line can be used to 
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constrain the neutron star parameters, let us see how it 
infers the black hole spin. The line shape, especially the 
extent of the red wing, carries the signature of the disc 
inner edge radius in the unit of stellar mass (n n c 2 /GM). 
This is because the red wing is primarily affected by lon- 
gitudinal Doppler effect due to the oribital motion of the 
disc matter, which broadens the line, and gravitational red- 
shift, which shifts the line towards lower energies. Since, 
for a black hole the disc can extend up to the innermost 
stable circular orbit (ISCO), the black hole spin parame- 
ter j — Jc/GM 2 (J: total angular momentum; M: mass), 
which determines the ISCO location, can be estimated by 
fitting the line shap e with an appropriate relativistic model 



fitting the line shap e with an appropriate relativistic model 
for Kerr spacetime (lLaor 199ll;lBeckwith and Don"e"ll2004l ; 



iDovciak et al. ll2004l : TMiller Il2007n " 

A neutron star system is usually more complex than 
a black hole syetem because of the following reasons. (I) 
While a spinning black hole is defined by only two parame- 
ters (mass and spin) , and the spacetime around it is the Kerr 
spacetime having analytical expressions, the structures and 
spacetimes of rapidly spinning neutron stars, usually har- 
boured by LMXBs, may have to be numerically calculated 
from at least two parameters apart from the EoS model 
l|Cook et al. Ill994h . (2) Unlike a black hole, a neutron star 
has a hard surface which can have observable effects. There- 
fore, for a neutron star, the disc m ay be terminated eithe r 
by ISCO or by the stellar surface l|Thampan et al.1 fl999). 
An example of the competition between these two e ffects 
has been shown in the Fig. 1 of Miller et al. I (|l998h and 
Fig. 1 of iBhattacharvva et all (|2000h . which demonstrate 
that usually circumferential radius n n of disc inner edge 
first decreases and then increases with the increase of stellar 
spin for a given EoS model and mass. This complicates the 
measurement of Jc/GM 2 and other stellar parameters us- 
ing the inferred circumferential n n c 2 /GM , if it is not known 
whether the disc terminates at ISCO or at the stellar surface. 
Therefore, theoretical computations of ri n as a function of 
stellar parameters for various EoS models are essential. Such 
computations will also be useful to determine if a measured 
r- ln c 2 /GM directly gives the neutron star radius-to-mass ra- 
tio Rc 2 /GM. Note that the computations of r- ln and stellar 
parameters will involve the numerical calculations of rapidly 
spinning neutron star structures, and hence Kerr spacetime 
cannot be used. Moreover, in order to constrain the EoS 
models, directly measur able neutron sta r para meters (e.g., 
Rc 2 /GM, M, i^spini S66 IBhattacharvva! (|201Ch ) should be 
expressed as functions of r ln c 2 /GM. Since Jc/GM 2 cannot 
usually be measured directly, a Jc/GM 2 vs. r ln c 2 /GM plot 
may not be very useful to constrain the EoS models. In this 
paper, such plots (Figs [1] [2] and [3} have been shown for 
comparisons with Kerr spacetime and to gain insight (§0). 

The procedure of constraining EoS models using di- 
rectly measurable parameters vs. r ln c 2 /GM relations can 
be utilized only if ri n c 2 /GM<6, because a much larger 
value of TinC 2 /GM might imply the truncation of the disc 
by other effects (see § |4}. Particularly useful would be 
rj n c 2 /GM < 6, because such values will confirm the effect of 
neutron star spin on a corotating disc. Therefore, a crucial 
question is whether observations show that r ln c 2 /GM<6. 
ICackett et al. I (|2010h report that r- ln (? /GM values fall into 
a small range of 6 — 15 for most of the neutron star LMXBs 
with established relativistic disc lines. Moreover, although 



these authors could not measure a value of r- ln c 2 /GM less 
than 6 (as they used non-spinning, i.e., Schwarzschild space- 
time), many of their fitted n n c 2 /GM values across the 
sources pegged at the lower limit 6. This happenned for both 
phenomenological and reflection models, and even when the 
Com pton broadening was taken into account l|Cackett et al. I 
2010). Such pegged best-fit values strongly suggest that 
n n c 2 /GM < 6 for many cases, which is supported by 
the disc line fitting with spinning Kerr spacetime models 
llBhattacharvva and Strohmayer ||2007l ; IPapitto et aL 1l2009l ; 
iReis et al. 1120091 ; iD'Af et al. Il2010h . This provides a good 
motivation to theoretically explore the above mentioned pro- 
cedure to constrain the EoS models. 

Since, in this paper, we have computed neutron star 
structures and the corresponding ISCO locations, let us 
briefly review some of t he pre vious studies on this topic. 
iKluzniak and Wagoner I l|l985l ) compu ted ISCO l o cation 
with slow stellar spin approximation. ICook et al. I (| 1994ft 
numerically calculated rapidly spinning neutron star struc- 
tures in full g enera l relativity using realistic EoS models. 
Miller et al. ( 1998) used these structure calculations to 
compu te ISCO locations for a few cases. IStergioulas et al~l 
( 1999) studied if strange stars in LMXBs could be ex- 
cluded using the orbital frequency at ISCO. The rele- 
vance of higher-order multip oles on the ISCO l ocatio n 
was analytically probed by IShibata and Sasaki I l|l998h . 
iBerti and Stergioulas I l|2004ft computed stellar models, mul- 
tipole moments and I SCO locations , and compared the 
ISCO results with the iManko et ai~l (|2000al lbh results for 
analytic solution for spacetime arou nd rapidly spinning 
neutro n stars, Kerr ISCO re sults and IShibata and~S asaki 
(|l998h results. The results of IBerti and Stergioulas I (|2004ft 
show that higher multipole moments cause differences be- 
tween rapidly s pinnin g neut ron star ISCO and Kerr ISCO . 
iPachon et all (120061 ) and ISanabria-Gomez et all (12010ft 
studied the ISCO a r ound a spinning magnetic neutron star. 
Ab ramowicz et al. I i|2003ft calculated circular geodesies in 
the Hartle-Thorne metric for slowly spinning neutron stars. 
IBerti et al. I (|2005ft found that Hartle-T horne approximation 
gives ISCO radii accurate to within 1%. Bej ger et al. I (|201dh 
suggested approximate analytic expressions for circular or- 
bits around rapidly spinning neutron stars. 

Disc lines have so far been used to measure 
the b lack hole spin (e.g., iBrenneman and Reynolds! 
([2006J)). It has also been proposed that these lines 
could be used to put an upper limit on the neu- 
tron star radius (IBhattacharvva and Strohmayer I 120071 ; 



Cac kett et al. |l2008h a nd to constrain other stellar param- 
eters ( Bhattac haryya I l2010ft . However, to the best of our 
knowledge, plausible methods to measure these parameters, 
and hence to constrain the EoS models in a systematic way 
using appropriate spacetimes have not been studied so far. 
The aim of this paper is to explore these methods. We have 
done so by calculating a huge number of neutron star struc- 
tures, and the corresponding ISCO locations, if such a lo- 
cation is outside the neutron star. Although, a large set of 
papers reported the studies on such structures and loca- 
tions (as indicated in the previous paragraph) , none of these 
aimed to probe neutron star parameters using iron lines. In 
§ § 13] and § 31 we describe our method, give the results 
and provide a discussion respectively. Note that, since the 
accretion disc is believed to be thin, we have not considered 
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nonequatorial orbits in our calculations. We have also not 
considered counterrotating orbits, because that would im- 
ply nscoc 2 /GM > 6, which could be easily confused with 
truncations caused by other effects (see §[4|. 



2 METHOD 

In this section, we briefly describe the procedure to compute 
rapidly spinning neutron star structures, and their equilib- 
rium sequences. The spacetime around such a star can be de- 
scribe d by the following m etric (using c — G = 1; iBardeen I 
l|l97(t l; ICook et al. I (|l994h ^: 



ds' 



-e 1+p dt z + e M (dr z + rW) + e 7_p r 2 sin' 6(d(j> - udtf , 



where the metric potentials 7, p, a, and the angular speed 
(w) of the stellar fluid relative to the local inertial frame 
are all functions of r and 9. For a given EoS model, 
and assumed values of stellar central density and polar- 
radius to equatorial-radius ratio, Einstein's field equations 
can be solved to find out r and 6 dependence of 7, p, 
a and uj, as well as to obtain the stable stellar structure 



dCook et al. Hl994l;lDatta et al. 1 


19981; iBhattacharvva et al. 1 


120001. 2001a bile ; Bhattacharvva 


2002h. This equilibrium so- 



lution can then be used to calculate bulk structure param- 
eters (e.g., M, R, J) of the spinning neutron star. Note 



that, henceforth, R (= 



j7e-Pe)/2 



) will denote equato- 



ri al circumferential ra dius of the neutron star (equation B6 
of ICook et al. I |l994|)). The equations of motion of a test 
particle in the spacet i me a round such a star are given in 
iThampan and Datta I (|l998T ). For example, the radial equa- 

(dr/dr) 2 =E 2 -V 2 , where, dr 



tion of motion is r 



is the proper time, E is the specific energy and a constant 
of motion, and V is the effective potential. V is given by 
V 2 = e 7+p [l + + 2ljEI - uj 2 1 2 , where I is the specific 

angular momentum, a constant of motion. We determine the 
radius of ISCO using the condition V, rr = 0, where a comma 
followed by one r represents a first-order partial derivativ e 
with respect to r and so on (|Thampan and Datta Ifl998h . 
Thus we can compute both circumferential radius of ISCO 
(denoted by nsco) and R, and set the theoretical value of n n 
to the larger of them (see §[TJ). For each EoS model, we have 
calculated v B pin sequences of equilibrium structures keeping 
fspin = constant, and changing other parameters. Since ^ sp i n 
is not an input for computing structures, a number of iter- 
ations are usually needed to compute the structure for a 
desired fspin value. The other sequences (e.g., M sequence, 
Rc 2 /GM sequence; figures of § [3} have been calculated by 
interpolations of v sp m sequences. 

We have used four representative EoS models of widely 
varying stiffness properties. This ensures sufficient gener- 
ality of our r esults. We briefly describe these models be- 
low. Model A (|Sahu et al. Ill993l ): This very stiff EoS model 
with maximum non-spinning mass M max m 2.59M0 is a field 
theoretical EoS for neutron-rich matter in be ta equilibrium 
based on the chiral sigma model. Model B (|Akmal et al. I 
QUI): This stiff EoS model with M max « 2.2OM is the Ar- 
gonne vis model of two-nucleon interaction, with the three- 
nucleon interaction (Urbana IX [UIX] model ) and the ef- 
fect o f relativistic boost corrections. Model C l|Baldo et al. I 
1 19971 ): This intermediate EoS model with M max « 1.79M 
is a microscopic EoS for asymmetric nuclear matter, derived 



from the Brueckner-Bethe-Goldston e manybody theory wit h 
explicit three-body terms. Model D (|Pandharipande Hl97ll ): 
This very soft EoS model with M max ~ 1.41Mq assumes an 
admixture of hyperons with the hyperonic potentials simi- 
lar to the nucleon-nucleon potentials, but altered suitably 
to represent the different isospin states. 



3 RESULTS 

We have computed v sp i n sequences (§ [2]) for 15 z^pm val- 
ues in the range of — 750 Hz for each EoS model 
About 16000 neutron star structures have been calcu- 
|:a1pd to establish our results, and we give example fig- 
ures in this section using a fraction of our computed num- 
bers. The code to compute these structures and nsco val- 
ues is well tested jDatta et al.l 19981; IThampan and Datta 



19981 ; IBhattacharvva et al. ll2000l . l2001al lbllc!: lBhattacharvva' 
2002). So far Kerr spacetime has been used to model the iron 
lines from spinning neutron star systems. Therefore, first we 
examine how our Jc/GM 2 vs. n n c 2 /GM plot deviates from 
the corresponding Kerr curve, in order to find out if Kerr 
calculations for iron lines can give acceptable constraints 
on neutron star parameters. Fig. [T] shows when the equa- 
torial radius of the neutron star is smaller than the ISCO 
radius, the deviation is relatively small, but is not negligi- 
ble, depending on the values of z^pm and M and the chosen 
EoS model. But when the neutron star equat orial radius is 
larger than the ISCO radius (see, for example. iMiller et al. I 
(1998)), the deviation is large because stellar equatorial ra- 
dius increases with the increase of v sp in, and such a situ- 
ation does not occur for black holes (Kerr spacetime). For 
example, for the EoS model A, the stellar equatorial ra- 
dius is greater than the ISCO radius for all z/ sp m values for 
M — IAMq, and hence the deviation is always very large. 
However, for M = 2.0Mq, the stellar equatorial radius is less 
than the ISCO radius for smaller values of u sp i n (Fig. [2]). In 
this case, the curves for both EoS models A and B are closer 
to the Kerr curve compared to these curves for Al = IAMq. 
However, for M = 2.OM0, stable neutron star structures do 
not exist for EoS models C and D. Therefore, from Fig.[T]and 
Fig. [2] we find that the deviation is more for (1) higher v spui , 
(2) lower M, and (3) stiffer Eo S models. The first two ef- 
fects can also be seen in Fig. 1 of IMiller et aTl j 19981 ). These 
effects are expected for n n c 2 /GM = Rc 2 /GM, because all 
these three points result in the increase of R. Let us now try 
to understand these points for n n c 2 /GM — rjscoc 2 /GM . 
The first point is understandable, because for z^ sp j n = 
the spacetime outside even a neutron star is Schwarzschild, 
which is the non-spinning special case of Kerr. The second 
point could be understood from the fact that for lower M 
the neutron star is usually less compact (that is the hard 
surface is farther from the centre), causing its spacetime to 
deviate more from that of a black hole. The third point may 
be explained from the lesser stellar compactness for a stiffer 
EoS model for given M and v sp in- However, we find that if 
M/M max is kept fixed instead of M, the differences among 
the Jc/GM 2 vs. r in c 2 /GM curves for various EoS models 
are small (Fig. [3j . This indicates that M max may be a suit- 
able parameter to characterize an EoS model. 

After finding that the Kerr spacetime is usually not 
good enough to model n n c 2 /GM, and since this spacetime 
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cannot be used for neutron star parameter calculation, wc 
have explored ways to constrain stellar EoS models using ap- 
propriate general relativistic computations for neutron stars 
(see §[2]). We have checked the relation between r ln c 2 /GM 
and a few stellar parameters, which can be measured from 
independent methods. These parameters are u Bp i n , Rc 2 /GM 
and M. Whenever v spm can b e measured (using regu lar pul- 
sations or burst oscillations; iBhattacharvva I (|201Ct )>. it is 
measured very accurately. Therefore, in the Figs. [4] [5] [6] 
and [7] we have fixed //spin • Rc 2 /GM and M, on the other 
hand, may be constrained in a range (usin g thermonuclear 
X-ray bursts, binary orbital motions, etc.; IBhattacharvva I 
i^Old )). and hence we have used them as dependent vari- 
ables in these figures. Figs. [4] and \5\ show Rc 2 /GM vs. 
n n c 2 /GM plots for two values of u ep in and four EoS mod- 
els. For a given v sp i n value, and if the stellar equatorial ra- 
dius is less than the ISCO radius, each EoS model traces 
a distinct curve. These curves, which are more separated 
from each other for higher v sp - m , can be used to constrain 
EoS models from the n n c 2 /GM value inferred from the iron 
line fitting, and the Rc 2 / GM value measured independently 
(see IBhattacharvva I |201(t 0. These figures show that even 
if Rc 2 /GM is not well constrained, a suitable upper limit 
of r in c 2 /GM can reject softer EoS models. If the stellar 
equatorial radius is larger than the ISCO radius, an oblique 
straight line is found for all EoS models, and an inferred 
rinC /GM value directly gives the Rc 2 /GM value. Figs. 0] 
and[5]clearly show the value of n n c 2 /GM (for a given u sp i n ), 
above which r- ln c? /GM can be used to directly infer the EoS- 
model-independent Rc 2 /GM value. Figs. [6] and [7] show even 
if, instead of Rc 2 /GM, M is kno wn from independent mea- 
surement ( Bhattac haryya 1 12010| ) ■ n n c 2 /GM inferred from 
iron line can be used to constrain the EoS models. Even for 
an unknown M, a suitable upper limit of r in c 2 /GM can be 
used to reject softer EoS models. Fig. [8] explores how for 
a known M value, n n c 2 /GA'/ inferred from iron line and 
independently constrained Rc 2 /GM can be used to con- 
strain EoS models. Similar result is shown in Fig. [9] where 
Rc 2 /GM is known and M is reasonably constrained. How- 
ever, for these two procedures (Figs.[8]and[9]), measurements 
of Rc 2 /GM and M appear to be more important than an 
inferred n n c 2 /GM . 



4 DISCUSSION 

In this paper, we explore ways to constrain neutron star EoS 
models by comparing the inferred values of n n c 2 /GM with 
the theoretical values. A n n c 2 /GM value may be inferred 
by f itting the relativistic disc lines with appropriate mod- 
els (|Bhattacharvva Il201ol ). We have shown that r in c 2 /GM 
calculated from Kerr spacetime is inadequate to distinguish 
between EoS models even when r lVL c? /GM = nscoc 2 /GM , 
and can largely differ from the correct value for rapidly spin- 
ning neutron star spacetime when n n c 2 /GM = Rc? /GM. 
We have numerically computed a huge number of n n c 2 /GM 
values, assuming that the disc terminates either at ISCO or 
at the stellar hard surface. This should be at least approxi- 
mately true for n n c 2 /GM< 6, although systematics may be 
introduced due to the effects of magnetic field and radiative 
pressure. Such systematics would imply that any inferred 
rinC 2 /GM value is basically an upper limit of nscoc 2 /GM 
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Neutron star angular momentum parameter 
vs. disc inner edge radius to stellar mass ratio 
( r inC 2 /GM). The solid curve is for Kerr spacetime (§0. Vari- 
ous symbols give the curves for different neutron star EoS models 
(§0 for M = 1.4Af©, and v Bpin ranging from Hz to 750 Hz. 
The negative slope of the curves implies that the disc terminates 
at ISCO, while the positive slope implies that it terminates at 
the stellar surface. This figure shows that the realistic r; n c 2 /GM 
values for neutron stars can significantly deviate from the Kerr 
values. 
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Figure 2. Neutron star angular momentum parameter 
(Jc/GM 2 ) vs. disc inner edge radius to stellar mass ratio 
(r in c 2 /GM). Similar to Fig. [TJ but for M = 2.0M Q . Note that 
stable neutron star structures do not exist for this high mass for 
the softer equation of state models C and D (§[2j. 



or Rc 2 /GM, which can still be used to reject softer EoS 
models (Figs. 2] [S] [5] and 0. These systematics can be re- 
duced by detailed modelling of these effects, as well as by 
independent observations. We have studied the relations be- 
tween rj n c 2 /GM and several directly measurable neutron 
star parameters ijBhattacharvva ll20irJ l. in order to estab- 
lish new ways to constrain EoS models. We have found that 
this iron line method could be effective to constrain EoS 
models, if the neutron star spin frequency is independently 
measured IBhattacharvva I l201Ch . This work is timely and 
important, as it provides motivation for future X-ray mis- 
sions, and because of the rapid progress in the disc line field 
via observations with XMM-Newton, Suzaku and Chandra. 
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Figure 3. Neutron star angular momentum parameter 
(Jc/GM 2 ) vs. disc inner edge radius to stellar mass ratio 
(r in c 2 /GM). Similar to Fig.[T] but for M = 0.95 X Af max . Here 
M max is the maximum mass for a non-spinning stable neutron 
star for a given EoS. This figure shows that a fixed M/M max 
gives similar curves, and hence M max may be a suitable parame- 
ter to characterize an EoS model. 
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Figure 4. Neutron 
(Rc 2 /GM) vs. disc 

(nnC 2 /GM) for various EoS models (§ for u S] 
Note that the oblique straight line portions of the curves are for 
r in c 2 /GM = Rc 2 /GM. This figure shows how a r in c 2 /GM value 
inferred from iron line can be used to constrain the EoS models 
for a known f ap i n value ( § [3j - 
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Figure 5. Neutron star equatorial-radius-to-mass ratio 
(R.c 2 /GM) vs. disc inner edge radius to stellar mass ratio 
( r inC 2 /GA/) for various EoS models (§0 for v ap i n = 600 Hz 
(similar to Fig. [4j . 





3 







2 


5 : 




2 


: 








a 


1 


5 : 




1 


: 







5 : 



!/ „, = 200 Hz 



/ A 




5.0 5.2 5.4 5.6 5.8 6.0 6.2 6.4 
r in c 2 /GM 

Figure 6. Neutron star mass (M) vs. disc inner edge radius to 
stellar mass ratio (r; n c 2 /GM) for various EoS models (§ [2]l for 
^spin = 200 Hz. Note that the oblique straight line portions of 
the curves are for r ln c 2 /GM = Rc 2 /GM. This figure shows how 
a n n c 2 /GM value inferred from iron line can be used to constrain 
the EoS models for a known v sp i n value (§|3J. 
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Figure 7. Neutron star mass (M) vs. disc inner edge radius to 
stellar mass ratio (r- ln c 2 /GM) for various EoS models (§ O for 
^spin — 600 Hz (similar to Fig. [6}. 
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Figure 8. Neutron star cquatorial-radius-to-mass ratio 
(Rc 2 /GM) vs. disc inner edge radius to stellar mass ratio 
{r in c 2 /GM) for various EoS models (§ |2j for M = 1.4M . 
Note that the oblique straight line portions of the curves are for 
ri n c 2 /GM = Rc 2 /GM. This figure shows, for a known M, how a 
r- ln c 2 /GM value inferred from iron line can be used to constrain 
the EoS models (§EJ). 
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Figure 9. Neutron star mass (M) vs. disc inner edge radius to 
stellar mass ratio (r; n c 2 /GM) for various EoS models (§ E| for 
Rc 2 /GM = 4.0. This figure shows, for a known Rc 2 /GM, how a 
ri n c 2 /GM value inferred from iron line can be used to constrain 
the EoS models (§|3J. 
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